Simultaneous dimerization and SU(4) symmetry breaking of 4-color fermions on the square lattice 
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Using infinite projected entangled pair states (iPEPS), exact diagonalization, and flavor-wave theory, we show 
that the SU(4) Heisenberg model undergoes a spontaneous dimerization on the square lattice, in contrast to its 
SU(2) and SU(3) counterparts, which develop Neel and three-sublattice stripe-like long-range order. Since the 
ground state of a dimer is not a singlet for SU(4) but a 6-dimensional irrep, this leaves the door open for further 
symmetry breaking. We provide evidence that, unlike in SU(4) ladders, where dimers pair up to form singlet 
plaquettes, here the SU(4) symmetry is additionally broken, leading to a gapless spectrum in spite of the broken 
translational symmetry. 
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The search for new quantum phases of matter, one of the 
main themes in the field of electronic materials with strong 
correlations, has become an important aspect of the physics of 
cold atoms too. A recent development concerns the possibility 
of loading multi-color fermions in optical lattices, which has 
reopened the investigation of the SU(iV) Heisenberg model 
that describes the permutation of iV-color objects on a lattice 
[1—4]. In condensed matter, these models require fine tuning 
of parameters. For instance, the SU(3) Heisenberg model is 
realized in spin-1 models with equal bilinear and biquadratic 
interactions, while in standard antiferromagnetic materials the 
biquadratic interaction is much smaller than the bilinear one. 
In the same spirit, the SU(4) model can be seen as a symmetric 
version of the spin-orbital Kugel-Khomskii model [5], but in 
actual materials the orbital-orbital interaction is not rotation- 
ally invariant [6]. By contrast, the simple quantum permuta- 
tion embodied by the SU(iV) Heisenberg model is a realistic 
starting point to describe the Mott phase of ./V-color fermions 
at filling 1 /N (one particle per site). 

In this paper, motivated by the conflicting results published 
so far, we concentrate on the case of 4-color fermions on 
the square lattice. In the Mott insulating phase with one 
fermion per site and large on-site repulsion, the relevant ef- 
fective model is defined by the Hamiltonian 

where the transposition operator P,y exchanges two colors on 
neighboring sites, and J > 0. 

An investigation of the model based on a variational and 
mean-field study suggested a plaquette ground state [7], with- 
out a long-range color-order. Exact diagonalizations of 8 and 
16 site clusters reported the presence of low lying singlets [9], 
and, inspired by the spontaneous singlet plaquette formation 
in the SU(4) ladder [10], suggested that these low-lying sin- 
glets might live in the subspace of singlet plaquette coverings, 
with possibly plaquette long-range order [8]. A different con- 
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FIG. 1 : (Color online) Sketch of dimer and color order realized in the 
ground state of the SU(4) Heisenberg model on the square lattice. 



elusion has been reached in Ref. 11, where a gapless spin liq- 
uid with nodal fermions has been shown to be the variational 
ground state for a class of projected fermionic wave functions. 
The presence of a long-range color order with wave-vector 
(2ir/N,2ir/N) for N = 2 (Neel order) and N = 3 (three- 
sublattice order with diagonal stripes [13]) hints to the possi- 
bility of extending such an ordering to N = 4. Lastly, a chiral 
spin liquid is predicted based on a particular large-N limit, but 
not below N=5 [12]. 

In this Letter, we show that none of these possibilities is 
realized, and that in fact both spatial and SU(4) symmetries 
are simultaneously broken: the ground state is spontaneously 
dimerized, and on each strong dimer, two colors dominate 
and combine into a 6-dimensional irrep of SU(4). In the re- 
sulting pattern, the same color does not appear on neighbor- 
ing dimers, so that if two colors are dominant on one dimer, 
nearest-neighbor dimers are dominated by the remaining two 
colors, resulting in the Neel-like arrangement sketched in 
Fig. 1. 

The first hint that simple color ordering is probably not real- 
ized appears at the level of linear flavor- wave theory (LFWT), 
the very method that confirmed the presence of stripe-order in 
the SU(3) case [13]. The starting point is a Hartree approx- 
imation based on a site-factorized variational wave function. 
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FIG. 2: (Color online) Sketch of some color and bond-energy con- 
figurations found by LFWT (a-c) and iPEPS (d-e). In all plots, the 
thickness of the lines is proportional to the square of (Pij) on the 
corresponding bond. The color content reflects that of the Hartree 
starting point for LFWT, and the actual result for an appropriate, in- 
finitesimal symmetry breaking field for iPEPS. (a) 'Plaquette' state 
selected by quantum fluctuations within LFWT. Its energy per site 
Eq/N = —3 J/2, (b) 'Dimerized' state with LFWT energy per 
site Eo/N = —1.293 J. (c) Stripe-ordered state, with much higher 
LFWT energy per site E /N = (4/tt — 2)J ~ -0.727J. (d) iPEPS 
results with D = 2 and a unit cell 4x4. At this order, the bond- 
energy pattern of the ground state is similar to that of the plaquette 
state selected by LFWT. (e) iPEPS results for D = 12 and a unit 
cell 4x2. The bond-energy pattern of this state is similar to that of 
the dimerized state of LFWT. Two types of dimers A and B with en- 
ergy much lower than vertical (V) and horizontal (H) dimers can be 
identified. They differ from each other by their local color densities. 



At this level, any state with different nearest-neighbor colors 
minimizes the energy, and, as in the case of SU(3), there is a 
macroscopic large number of such states. The next level of ap- 
proximation consists in including quantum fluctuations within 
LFWT: The classical energy of each state is corrected by the 
zero point energy of a quadratic bosonic Hamiltonian which, 
on a nearest-neighbor bond with color a on site i and color b 
on site j, is given by "H fw = A^A^ - 1, with = a] + 6 4 
where a,j and 6j are bosonic operators. Since there is an in- 
finite number of Hartree ground states, we have used the fol- 
lowing strategy: compare simple, periodic states, and make 
a systematic exploration of all states on finite clusters. The 
resulting picture is summarized in Fig. 2(a-c). First of all, 
and most importantly, the state with stripe order of Fig. 2(c), 
which breaks SU(4) symmetry but leaves all nearest neighbor 
bonds equivalent, is clearly not favored. In fact, many peri- 
odic states which break the spatial symmetry between bonds 
concomitantly with the SU(4) symmetry have much lower en- 
ergy, the lowest one being a quadrumerized state in which 
strong bonds build 4-site plaquettes (Fig. 2(a)). Among pos- 
sible states we can also find a dimerized state (Fig. 2(b)) with 
energy slightly above that of the plaquette state but far below 
that of the stripe state. 

This result, rather surprising in view of the SU(3) case [13], 
can be rationalized as follows: i) In a given ground state, the 



LFWT Hamiltonian is the sum of independent Hamiltonians 
that describe the motion of bosons on the connected clusters 
spanned by pairs of colors. For example, if a pair of neigh- 
boring sites with say colors a and b is such that all nearest 
neighbor sites have colors different from a and b, it forms a 2- 
site connected cluster. For the state of Fig. 2(a), there are two 
types of clusters: 2-site clusters on strong bonds, and 4-site 
clusters (plaquettes) on weak bonds; ii) The zero point energy 
per bond tends to increase with the cluster size. In particular, 
on a 2-site cluster, the ground state energy of the Hamiltonian 
is equal to —1, i.e. there is no zero-point contribution to the 
energy [14], whereas larger clusters lead to finite frequencies, 
hence to strictly positive contributions to the zero-point en- 
ergy. In the stripe-ordered state of Fig. 2(c), all clusters are 
infinite, one dimensional zig-zag stripes, and no bond is able 
to minimize the energy. By contrast, half the bonds have en- 
ergy — 1 in the state of Fig. 2(a) (thick bonds) and a quarter of 
them in the dimer state of Fig. 2(b). The difference betweeen 
the present case and SU(3) comes from the fact that it is im- 
possible to realize isolated two-site clusters with only three 
colors on the square lattice, hence to have local bonds with 
very low energy. 

So flavor-wave theory gives strong indication that the 
ground state is not a simple color-ordered state, and that some 
kind of additional lattice symmetry breaking takes place. One 
may be tempted to go one step further, and to conclude that 
it predicts plaquette order. However, if the symmetry is bro- 
ken and small clusters appear, LFWT becomes inaccurate (for 
instance, it cannot restore the SU(4) symmetry in the ground 
states of a 4-site cluster), and one should rather turn to alter- 
native methods. 

Infinite projected entangled-pair states (iPEPS) results — 
Next, we present results obtained with iPEPS [15], general- 
ized to arbitrary unit cells [16], a variational approach based 
on a tensor network ansatz aimed at efficiently representing 
the ground state of two-dimensional lattice models in the ther- 
modynamic limit. It consists of a network of rank-5 tensors, 
one for each lattice site, where each tensor is connected to its 
nearest neighbors by four bond indices with a certain bond di- 
mension D, and the fifth index carries the local Hilbert space 
of a lattice site. The accuracy of the ansatz can be system- 
atically controlled by the bond-dimension D. We use Nt 
different tensors, arranged in a rectangular unit cell of size 
L x x L y — Nt which is periodically repeated in the lattice 
[16]. The tensors are optimized by performing an imaginary 
time evolution, using the so-called simple update [17, 18]. 
For the computation of expectation values we use the corner- 
transfer matrix method [16, 20] to approximately contract the 
tensor network, where the accuracy of the contraction can be 
controlled by another parameter called boundary dimension \- 
The simulation results presented in the following are extrap- 
olated in x, where the extrapolation error is small compared 
to the symbol sizes. To improve the efficiency for large bond 
dimensions we use tensors with Z g symmetry [21] [a discrete 
subgroup of SU(4)]. 

Figure 3(a) shows the iPEPS results for the energy per site 
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FIG. 3: (Color online) iPEPS results as a function of the inverse bond 
dimension 1/D. Dotted lines are only a guide to the eye. (a) Energy 
per site in units of J for different unit cell sizes in the iPEPS, com- 
pared to the variational Monte Carlo (VMC) energy from Ref. [11] 
(extrapolated to the thermodynamic limit), and exact diagonalization 
results, (b) Color order parameter (COP), which is zero if on each 
site the color density of each color k is the same, i.e. = 0.25. 
The data does not seem to extrapolate to zero in the limit D — > oo 
and we therefore expect the SU(4) symmetry to be spontaneously 
broken, (c) Bond energies of the bonds H, V, A and B marked in 
Fig. 2 obtained from the 4x2 unit cell for D > 4 and from the 4x4 
unit cell for D < 4. The dimerization (strong A/B bonds) is seen to 
persist in the limit D — > oo. 



as a function of 1/D for different unit cell sizes. The energy 
obtained with the 2x2 unit cell is considerably higher than 
the one obtained with the 4x2 and 4x4 unit cells. For D < 4 
the lowest energy is obtained with the 4x4 unit cell, while 
for D > 4 similar results are obtained with the 4x2 unit 
cell. This indicates that the ground state breaks translational 
invariance in a way which is compatible with the 4x2 unit 
cell, but not with the 2x2 unit cell. For a bond dimension 
D — 8 the iPEPS variational energy is comparable to the vari- 
ational Monte Carlo (VMC) energy from Ref. [11], however, 
for larger D we obtain energies which are considerably lower. 
We note that the energy has not converged yet with bond di- 
mension D, and thus further corrections to the ground state 
energy (and other quantities) can be expected for larger D. 

We next turn our focus to local properties of the ground 
state. The bottom part of Fig. 2 shows the local color den- 
sities on each site and the energy on each bond for D = 2 
(left) and D — 12 (right). The evolution with D of the en- 
ergy of the various bonds is shown in Fig. 3(c). For D = 2, 
the bond pattern is similar to that of the linear FWT ground- 
state of Fig. 2(a). However, as soon as D > 2, three types of 
bonds appear in the ground state, and the bonds form a colum- 
nar dimerized pattern similar to the linear FWT of Fig. 2(b). 
Since increasing D allows more quantum fluctuations, we in- 
terpret this result as an indication that the dimerized pattern 
is ultimately stabilized when enough quantum fluctuations are 
taken into account. 

For SU(2) antiferromagnets, if spontaneous formation of 



dimers occurs, this is the end of the story: The ground state 
degeneracy is given by the number of equivalent dimer cov- 
erings, and each ground state is adiabatically connected to a 
product of singlets on the dimers. However, for SU(4), it takes 
at least 4 sites to make a singlet, and on a dimer, the ground 
state is 6-fold degenerate. This 6-dimensional subspace cor- 
responds to the 6-dimensional irrep of SU(4). A convenient 
basis is provided by the states (\a(3) — |/?a))/v2 where the 
couple (a, j3) is a pair of different colors chosen out of the col- 
ors of the model and can take 6 values. These states are con- 
nected by SU(4) rotations. So, even in the presence of sponta- 
neous dimerization, there is room for further SU(4) symmetry 
breaking. In our iPEPS results for D = 12, one can immedi- 
ately recognize the formation of two types of dimers, marked 
with A and B in Fig. 2(e), which are arranged in a (colum- 
nar) checkerboard order. The energy of both A and B dimers 
is identical, however, the local color densities are different. 
In dimers A the first two colors [blue and yellow in Fig. 2(e)] 
are dominant, whereas in the B-dimers the density of the other 
two colors is larger. [Here we applied a small initial field to fix 
the "direction" of the SU(4) symmetry breaking in the space 
of the four colors. We note that the same results (up to SU(4) 
rotations) are obtained without an initial field]. This type of 
order is seen to persist also in the limit D — > oo, as shown 
in Fig. 3(b). We also verified that in the simulations done 
with a 4 x 4 unit cell this pattern is repeated twice. Thus, 
our iPEPS results suggest that both the SU(4) symmetry and 
translational symmetry are broken, resulting in a Neel-like or- 
der with dimers alternating between pairs of colors. 

Exact diagonalization (ED) results — We now compare 
these predictions with the results of ED of finite clusters. The 
energies per site \E/(NJ)] in the ground state are —1 for 
N = 8, -1.0844 for N = 16 [9] and -1.0008 for N = 20. 
In Fig. 4 we present results for several correlation functions 
obtained for the N = 16 system. In Fig. 4(a) Fourier- 
transformed color-density correlations are shown. The diam- 
eters of the circles are proportional to the Fourier component. 
The maxima are found at momentum (it, tt/2) and equivalent 
positions. This provides a possible explanation as to why the 
N = 8 and N = 20 samples have higher energy per site 
than the N = 16 cluster, as their Brillouin zone does not con- 
tain these particular momenta and thus frustrate the preferred 
order [22]. In Fig. 4(b) we display the connected nearest- 
neighbor bond energy correlations (PijPki) — (Pij) 2 - The 
correlations present a clear striped signal and are in qualitative 
agreement with the bond energy modulation pattern obtained 
using iPEPS, however the modulation is not very strong. Fi- 
nally in Fig. 4(c) we calculate correlations between the pro- 
jectors onto a given basis state of the d = 6 representation 
of SU(4) on a bond (Vta^ih j)'P(a,V)(f : > 0)- This correlation 
supposedly captures the Neel structure on top of a columnar 
dimerized background. Indeed the correlations in (c) show the 
expected four positively correlated bonds at the correct loca- 
tion confirming the qualitative picture put forward in Fig. 1. 
We have also measured the spin chirality correlations advo- 
cated in Ref. [12]. We indeed find correlations pointing to- 
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FIG. 4: (Color online) (a) Color structure factor for the N = 16 square sample in momentum space. The maxima of the structure factor (filled 
circles) are located at momentum (it, n/2) and equivalent points, (b) Connected bond energy correlations (P lJ P k i) - (Pij) 2 for N = 16. 
(c) (P( a ,b){i, j)P(a,b)(k, I))- In (b) and (c) the line width is proportional to the absolute value of the correlations and blue, full (red, dashed) 
lines denote positive (negative) values, (d) Low energy spectrum of the N = 16 sample plotted as a function of the quadratic SU(4) Casimir 
operator C. Expected levels in the tower of states of the SU(4) symmetry breaking state sketched in Fig. 1 are highlighted by dashed circles. 
In the grey shaded region there are additional levels (not shown) which have not been resolved as function of the Casimir. 



wards uniform chirality on triangles, however the correlations 
decay quickly over the size of the sample, implying the prob- 
able absence of a chiral spin liquid phase on the SU(4) square 
lattice. 

We conclude this ED section by presenting the low energy 
tower of state in Fig. 4(d) as a function of the quadratic SU(4) 
Casimir operator. In the case of either discrete or continu- 
ous symmetry breaking the low energy spectrum is expected 
to exhibit characteristic energy levels which become degen- 
erate with the ground state in the thermodynamic limit and 
which enable the spontaneous symmetry breaking. For the 
proposed state in Fig. 1 we expect a tower of states which 
is a fourfold superposition of the tower of state of a SU(4) 6 
bipartite Neel state. The later is given by the irrep content 
C = 0, 4, 6, 10, 12, . . .. Combining with the spatial degener- 
acy associated with the four distinct dimer configurations, the 
even members of the Casimir sequence are expected to be- 
long to the spatial symmetry sectors (0, 0)Ai, (0, 0)i?i and 
the twofold degenerate sector (tt, 0)a x = — 1, a y = 1, while 
the odd members belong to the fourfold degenerate momen- 
tum (7T,7r/2) a x = 1 sector. When comparing this predic- 
tion (large circles) to the actual low energy spectrum shown 
in Fig. 4(d) one observes that the lowest levels for C — and 
4 are indeed the ones expected for the proposed state in Fig. 1 . 

Conclusion — The various approaches followed in this 
paper point to a rather coherent picture: a spontaneous dimer- 
ization appears in the ground state, and the remaining degrees 
of freedom, an irrep of dimension 6 on each dimer, develop 
Neel-like order. The model that describes the fluctuations 
around one of the dimerized ground states should thus be a 
model with SU(4) symmetry with the 6-dimensional irrep at 
each site of a square lattice and with anisotropic couplings - 
horizontal and vertical pairs of nearest-neighbor dimers are 
inequivalent. In that respect, it is interesting to note that the 
isotropic version of the SU(4) Heisenberg model with the 6- 
dimensional irrep and only nearest-neighbor couplings has 
been studied by Quantum Monte Carlo [23] and variational 



Monte Carlo [24] with conflicting conclusions. While in the 
latter study long-range Neel order is found, the former pre- 
dicts that the correlations decay algebraically, i.e. that there is 
no true long-range order. Coming back to the evidence in fa- 
vor of Neel order in our case, the presence of short-range cor- 
relations is clear in view of the ED and iPEPS results. How- 
ever, correlations revealed by ED are not very strong, and the 
prediction of long-range order relies mostly on the empirical 
1 /D scaling of the iPEPS order parameter of Fig. 3. It would 
thus be very interesting to challenge this prediction by deriv- 
ing and studying the effective model relevant to the present 
case. 
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